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Summary 

The National Aeronautics and Space Adminis- 
tration, at its Lewis Research Center, develops and 
coordinates plans and results related to icing 
computer-code predictions, icing research tunnel 
(IRT) measurements, and icing research flights. To 
allow for the comparison of research flight results 
with icing code predictions and IRT measurements, 
Lewis has collaborated with the NASA Langley 
Research Center for the planning and analysis of a 
series of icing research flights related to the deter- 
mination of stability and control derivatives from re- 
search flight data. 

This paper presents the results of applying a mod- 
ified stepwise regression algorithm and a maximum 
likelihood algorithm to flight data from a twin-engine 
commuter-class icing research aircraft. The results 
are in the form of body-axis stability and control 
derivatives related to the short-period, longitudinal 
motion of the aircraft. Data were analyzed for the 
baseline aircraft ( “uniced” ) and for the aircraft with 
an artificial glaze-ice shape attached to the leading 
edge of the horizontal tail. The results are discussed 
as to the accuracy of the derivative estimates and 
the difference between the derivative values found for 
the baseline and the “iced” aircraft. Additional com- 
parisons are made between the maximum likelihood 
results and the modified stepwise regression results 
with causes for any discrepancies postulated. 

Introduction 

The known dangers to safe flight that are caused 
by the accumulation of ice on aircraft components 
have given rise to a multifaceted approach to un- 
derstanding the causes and predicting the effects of 
aircraft icing (ref. 1). The National Aeronautics 
and Space Administration (NASA), at its Lewis Re- 
search Center (LeRC), supports an Icing Technology 
Project Office in which plans and results related to 
icing computer-code predictions, icing research tun- 
nel (IRT) measurements, and icing research flights 
are coordinated. To allow for the comparison of re- 
search flight results with icing code predictions and 
IRT measurements, LeRC has collaborated with the 
NASA Langley Research Center (LaRC) for the plan- 
ning and analysis of a series of icing research flights 
related to the determination of stability and control 
derivatives from research flight data. 

The successful estimation of stability and control 
derivatives from flight data has two important ram- 
ifications. First, the values for the derivatives can 
be compared with values derived from the analytical 
icing codes and from the IRT. These comparisons al- 
low for an assessment of the confidence that should be 


put in analytical predictions and wind-tunnel results 
as they relate to an aircraft in flight. Second, flight- 
derived derivatives can be used judiciously along with 
those provided by analytical predictions and wind- 
tunnel tests to upgrade simulator math models to 
provide a realistic set of aerodynamics for pilot-in- 
the-loop simulations of icing scenarios. Such simula- 
tions can provide definitive information on how any 
measured degradation of the aircraft flying qualities 
actually affects its handling qualities. 

Previous work on the estimation of stability and 
control derivatives from flight tests for an iced air- 
craft has been performed by LeRC and Kohlman 
Systems Research, Inc. (KSR), as reported in refer- 
ences 2 and 3. That earlier work employed a modified 
maximum likelihood technique to estimate the stabil- 
ity and control derivative values. The accuracies of 
these stability and control derivatives were estimated 
by an approximation to the Cramer-Rao bound given 
by the maximum likelihood algorithm. Since it is 
known that the Cramer-Rao bound estimates indi- 
cate better accuracy than is actually achievable by 
repeated experiments (ref. 4), and since the earlier 
work showed small, but apparently discernible, dif- 
ferences in the derivative values for the “iced” versus 
“uniced” aircraft, the first phase of the LaRC plan 
was to determine the expected accuracy of the esti- 
mated derivative values. To this end an ensemble of 
45 maneuvers was flown at the same flight condition. 
These maneuvers were then individually analyzed as 
to longitudinal stability and control derivatives us- 
ing both a maximum likelihood (ML) algorithm and 
a stepwise regression (SR) algorithm. The ensem- 
ble standard deviation for the parameter estimates 
was compared with the average of the Cramer-Rao 
bound estimates (maximum likelihood) for the en- 
semble and the average estimate of standard error 
(stepwise regression). Maneuvers were then flown 
from constant-power, 1 g flight at several trim air- 
speeds for both the uniced and artificially iced air- 
craft, and the recorded data from those maneuvers 
were analyzed for longitudinal stability and control 
derivatives using both a stepwise regression and a 
maximum likelihood algorithm. The flight test pro- 
gram also included several data compatibility ma- 
neuvers for the assessment of data quality and several 
deceleration/acceleration maneuvers for direct deter- 
mination of lift and pitching- moment curves. 

The purpose of this paper is to present the re- 
sults of applying a modified stepwise regression algo- 
rithm (ref. 5) and a maximum likelihood algorithm 
(ref. 6) to flight data from a twin-engine, commuter- 
class icing research aircraft. The results are in the 
form of body-axis stability and control derivatives 
related to the short-period, longitudinal motion of 
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the aircraft. Data were analyzed for the baseline 
aircraft (“uniced”) and for the aircraft with an ar- 
tificial glaze-ice shape attached to the leading edge 
of the horizontal tail. The results are discussed as 
to the accuracy of the derivative estimates and the 
difference between the derivative values found for the 
baseline “uniced” and the “iced” aircraft. Additional 
comparisons are made between the maximum like- 
lihood results and the modified stepwise regression 
results with causes for any discrepancies postulated. 
The paper is organized as follows: After this intro- 
duction, the aircraft, its instrumentation, and the 
flight test maneuvers will be briefly discussed fol- 
lowed by a section on the analysis techniques. The 
main part of the paper will then present the results of 
the data analysis followed by a section summarizing 
the conclusions to be drawn from these results. 
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Symbols 

a z vertical acceleration, g units 


s e{ ) 


cost function for maximum 
likelihood algorithm 

cost function for stepwise 
regression algorithm 

aircraft mass, kg 

number of data points 

roll, pitch, and yaw rates, re- 
spectively, rad/sec or deg/sec 

wing area, m 2 

total sum of squares (corrected 
for mean) 

standard error 

average of standard errors 
given by program 

ensemble-average standard 
error for a group of values 


b 

bax 

baz 

b q 


be 


Ci,C m ,C n 


Cx^Y,C Z 


c 

9 


IxJyJz 


txz 


wing span, m 

calculated constant bias in 
longitudinal- acceleration 
measurement, m/sec 2 

calculated constant bias in 
vertical-acceleration measure- 
ment, m/sec 2 

calculated constant bias in 
pitch-rate measurement, 
rad/sec 

calculated constant bias in 
angle-of-attack measurement, 
rad 

calculated constant bias in 
pitch-angle measurement, rad 

rolling-, pitching-, and 
yawing- moment coefficients, 
respectively 

longitudinal-, lateral-, and 
vertical-force coefficients, 
respectively 

mean aerodynamic chord, m 

acceleration due to gravity 
(1# « 9.81 m/sec 2 ) 

moments of inertia about 
roll, pitch, and yaw axes, 
respectively, kg-m 2 

product of inertia, kg-m 2 


^«*/2,n— 2 

student ^-distribution 

u, v , w 

longitudinal, lateral, and 
vertical velocity components, 
respectively, m/sec 

V 

total airspeed, (n 2 -b'c 2 +'ie 2 ) 1 / 2 , 
m/sec 

xj(i) 

measured aircraft response or 
control surface input 

y(i) 

computed aerodynamic 
coefficients 

yo( x o) 

linear approximation of Cm 6 
corresponding to some xq 
value 

z i 

measured output 

h 

computed output 

a 

corrected (free-stream) angle 
of attack, rad or deg 

a M 

experimentally measured angle 
of attack, rad or deg 

a* 

criteria for achieving a desired 
prediction interval 

0 

angle of sideslip, rad or deg 

bei^a,) b r 

elevator, aileron, and rudder 
deflection, respectively, rad or 
deg 

0 

aerodynamic stability or 
control parameters 
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00 


constant offset term 

©i 


derivative with respect to ith 
response or control surface 
input 

9,<f> 


pitch and roll angles, respec- 
tively, rad or deg 

A a 


upwash (scale-factor) correc- 
tion for angle of attack , 

P 


air density, kg/m 3 
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variance estimate of measure- 
ment noise (eq. (2)) 
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A dot over a symbol represents the derivative with 
respect to time. 

Aircraft and Instrumentation 


the 24-sec slot while retrimming and/or housekeeping 


Research personnel onboard the aircraft were the 


in the cabin for operating the DAS and coordinating 
available 24- sec test slots with the pilot. 

After each flight, the recorded data cassettes were 
t to KSR where the data were reduced to engineer- 
ing units. The engineering-units data were written by 
KSR to a nine-track magnetic tape which was sent 


Research Flights 

The four airplane configurations flown are given 


Configuration 

Aircraft condition 

Flap 

deflection 

Flight 

designation 

1 

Baseline (uniced) 

0° 

88-5 

2 

Baseline (uniced) 

10° 

88-6 

3 

Artificial ice on tail 

0° 

88-7 

4 

Artificial ice on tail 

10° 

88-8 


The icing research aircraft is a modified deHavil- 
land DHC-6 Twin Otter. It is powered by two Pratt 
& Whitney PT6A-20A gas turbine engines. Physical 
characteristics for the aircraft are found in table I. A 
three- view drawing is presented in figure 1. Longi- 
tudinal control is effected by elevator and hydrauli- 
cally actuated flaps; lateral and directional control, 
through ailerons and rudder, respectively. The ele- 
vator, ailerons, and rudder all have trim tabs. 

The instrumentation, which was provided and in- 
stalled by KSR, included three- axis orthogonal lin- 
ear accelerometers, three-axis angular rate gyros, 
and three-axis attitude gyros. Angles of attack and 
sideslip were measured by balsa wood vanes installed 
on a nose boom. Elevator, rudder, and aileron posi- 
tions were measured by potentiometers placed at the 
control surface. 

Since flaps were set and remained constant during 
an entire series of maneuvers, flap position was given 
by a readout in the cockpit and was recorded by the 
flight engineer for future reference. 

Data were recorded onboard at 20 samples per 
second; all channels were sampled within 1 msec by 
the KSR data acquisition system (DAS). Because of 
buffer limitations of the KSR system, only 24 sec 
of the 20-samples-per-second data could be obtained 
before the buffer contents had to be written onto a 
cassette tape. The write-off period was nominally 60 
to 90 sec. Hence, maneuvers had to be flown within 


These flights were designated as 88-5, 88-6, 88-7, 
and 88-8, respectively, by LeRC. Within each of these 
flights, each 24-sec slice of individual maneuvers was 
assigned a run number. The aircraft is limited to 10° 
flaps in an icing environment as a safety precaution. 
The artificial ice was a styrofoam form modeled 
after moderate glaze ice that had previously been 
seen and carefully photographed during natural icing 
encounters and in the IRT. This shape was glued 
to an aluminum form that was contoured to the 
leading edge of the horizontal tail. The aluminum 
was secured to the leading edge by attaching straps 
to the stabilizer hinge line as shown in figure 2. Four 
types of maneuvers were flown in each configuration: 

1. Data compatibility maneuvers 

2. Small-amplitude longitudinal perturbations 

3. Large-amplitude longitudinal perturbations 

4. Deceleration/acceleration maneuvers 

The data compatibility maneuvers involved excit- 
ing large angular rates about all axes, thus leading 
to large excursions in all Euler angles and to large 
variations in linear accelerations, airspeed, angle of 
attack, and angle of sideslip. These data were ex- 
amined using the algorithm described in reference 7 
to assess the quality of the DAS measurements. In 
these maneuvers, in which no attempt was made to 
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excite any of the aircraft modes, only the air data 
sensors, rate sensors, and linear accelerometers were 
exercised. 

The small-amplitude maneuvers were initiated by 
elevator “2-1” doublets designed to excite the short- 
period longitudinal motion of the aircraft. An ex- 
ample of the control surface movement and aircraft 
longitudinal responses is given in figure 3. These 
maneuvers were planned to excite ±0.3# of verti- 
cal acceleration. However, because of an erroneous 
accelerometer readout in the cockpit, these maneu- 
vers were sometimes of amplitude ±1.0# of vertical 
acceleration. 

The large-amplitude perturbations were initiated 
by random elevator doublets superimposed on an 
otherwise monotonically increasing angle of attack. 
For these maneuvers, the pilot pulls steadily back on 
the stick over the 24-sec test window, interrupting 
the steady pull at random intervals to apply small 
elevator doublets. In this way a large angle-of- 
attack range can be probed in a single maneuver at 
a constant power setting. 

The deceleration/acceleration maneuvers were 
performed by setting power for 1.0^-level flight at 
120 knots and then slowly pulling back on the stick, 
decelerating at about 2 knots/sec. Then, if time 
was available in the 24-sec test window, the pilot 
would push the stick forward again to accelerate at 
2 knots/sec, eventually returning to the initial flight 
condition. 

The small-amplitude maneuvers fell into two 
groups. In the first group, to ascertain the accu- 
racy of the derivative estimates, a set of maneuvers 
was performed in which the aircraft was trimmed at 
an airspeed of approximately 120 knots, and a series 
of 45 small-amplitude longitudinal perturbations was 
executed. These 45 sets of small-amplitude maneu- 
vers provided a statistically significant ensemble of 
maneuvers, each of which should produce the same 
derivative estimates. Analyzing each maneuver of 
this ensemble and calculating the mean and standard 
deviation for the 45 results for each of the longitudi- 
nal stability and control derivatives gives a measure 
of how much a parameter estimate might vary sta- 
tistically between two maneuvers at the same flight 
condition. 

In the second group, the rest of the small- 
amplitude maneuvers were performed at different air- 
speeds with one or two repeats at each airspeed. 
These maneuvers were analyzed for trends in the 
derivatives relative to flight condition. Because the 
effect of artificial ice on lift and drag of the whole 
aircraft should be small and because the aircraft was 
trimmed at different airspeeds using only the eleva- 
tor tab (constant power), it was assumed that direct 


comparisons could be made between the clean and 
“iced” derivatives at each trim airspeed. 

Data Analysis Methods 

The raw engineering- units data are first corrected 
for all known instrumentation offsets from the air- 
craft center of gravity and for flow effects. Next, 
several maneuvers are analyzed with the maximum 
likelihood algorithm reported in reference 7 for the 
estimation of instrument biases and possible scale- 
factor errors. This process is known as data compat- 
ibility or data consistency checking. Once a data set 
is available, the identification process still has two 
remaining stages. The first stage is the determina- 
tion of the parametric model structure. Then, the 
second stage is the estimation of the corresponding 
parameter values. The model structure was deter- 
mined using a modified stepwise regression. This 
technique can determine significant terms among 
the candidate variables and estimate corresponding 
parameters. Each new variable chosen for entry into 
the regression equation is the one that has the largest 
correlation with the measured dependent variable 
(y) after adjusting for the effect on y of the vari- 
ables already selected. This correlation is reflected 
by the calculated F-statistic of the variable. The se- 
lected parameters are estimated by minimizing the 
cost function 



where y(i) is the aerodynamic coefficient computed 
from measured data at time i, N is the number of 
data points, and t + 1 is the number of parameters 
(or independent variables) in the regression equation. 
At every step of the regression, the variables incorpo- 
rated into the model in previous stages and the new 
variable entering the model are reexamined. Any 
variable that provides a nonsignificant contribution 
(due to correlation with more recently added terms) 
is removed from the model. The process of select- 
ing and checking variables continues until no more 
are admitted to the equation and no more are re- 
jected. Experience shows, however, that the model 
based only on the significance of individual parame- 
ters in the model equation can still include too many 
terms and, therefore, may have poor prediction ca- 
pabilities (ref. 8). Therefore, several criteria for the 
selection of an adequate model were introduced in 
reference 8. Details of this model-structure determi- 
nation procedure are also explained in that reference. 

Finally, each model and each maneuver were an- 
alyzed by the maximum likelihood (ML) method 
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described in reference 6 assuming no external noise. 
The corresponding cost function has the form 

N k 

JmL = ^2 \ z j^) ~ %(*’ ©)] ( 2 ) 

i= i i=i 

where k is the number of output variables, is 
the variance estimate of measurement noise, zj and 
Zj are the measured and computed outputs, respec- 
tively, and © is a vector of estimates of unknown 
parameters. 

In addition to the parameter estimates, the ML 
method provides an estimate of the accuracy of each 
parameter through an approximation to the Cramer- 
Rao bound for that parameter. Theoretically, the 
Cramer-Rao bound presents a lower bound for the 
variance of the estimates of a parameter. The ap- 
proximation to the Cramer-Rao bound is realized by 
the diagonal elements of the inverse of the Fisher in- 
formation matrix. Since with flight data there is no 
“true” value with which to compare parameter esti- 
mates, agreement between the estimates of a given 
parameter from the ML and SR algorithms increases 
the overall confidence in that estimated value. 

Results and Discussion 

Data Compatibility 

Three maneuvers that were designed specifically 
for a longitudinal data-compatibility check were 
analyzed using the algorithm developed in refer- 
ence 7. This compares measured and predicted air- 
craft responses using the kinematic equations as a 
model. It employs the maximum likelihood tech- 
nique for finding uncorrected systematic errors in 
measured responses in terms of constant bias errors 
and scale-factor errors. For the symmetric longi- 
tudinal maneuvers, possible biases in angle of 
attack (b a ) y pitch rate (b q ), pitch angle (&#), lon- 
gitudinal acceleration (b ax ), and vertical accelera- 
tion ( b az ) were estimated. The resulting values 
are found in table II. In addition to the six bias 
values estimated, a scale factor in angle of at- 
tack (A a ) was determined as a correction for up- 
wash at the tip of the nose boom. Although 
computed angle-of-attack time histories appeared to 
correspond well with those measured in flight, the 
differences between measured and computed data 
(residuals) still exhibited a definite structure (were 
not randomly distributed). (See fig. 4.) As a possi- 
ble remedy for this situation, several candidate time 
shifts were introduced in the angle-of-attack time 
histories. It was postulated that this action might 


account for an uncorrected lag present in the data 
recording technique or in effects of the dynamics 
of the vane. Although increasing time shifts did 
“whiten” residuals, there was no physical justifica- 
tion for corrections of the magnitude of shift required. 

Because of these problems involving angle-of- 
attack values and an overall inconsistency in com- 
puted bias values for the three compatibility runs, the 
45 repeated longitudinal maneuvers and 3 additional 
large- amplitude maneuvers were analyzed using the 
algorithm mentioned above. From the analysis of 
these maneuvers, mean values of upwash correction 
(A a ) and all bias terms were compiled (table III). 

The magnitude and consistency of bias values ob- 
tained from analyzing all runs gave sufficient confi- 
dence in the data acquisition system for the small 
perturbation maneuvers exciting the short-period 
modes of the aircraft. They were not used as cor- 
rections to the data because of their very low values 
and relatively high standard errors, a result which 
indicates that no real biases exist. 

The average scale- factor estimate for upwash cor- 
rection from the 45 repeat maneuvers was 

A q = 0.10605 

with the Cramer-Rao bound estimate of 
3(A a ) = 0.00660 

and an ensemble standard deviation of 

s E (X a ) = 0.00741 

The upwash correction was made to all measured a 
in the following manner: 

a M = (1 + A a )a 

or 

a = (1 + A a )~ l a M 

where ajvf ls the experimentally measured angle of 
attack and a is the corrected free-stream angle of 
attack. 

Accuracy of Derivative Estimates 

The purpose of executing the 45 maneuvers under 
identical flight conditions was to acquire a series of 
flight data sets that would yield a statistically sig- 
nificant base of parameter and error estimates. The 
benefits of such an ensemble are twofold. First, the 
standard error calculated for the ensemble of param- 
eter estimates themselves gives us a measure of the 
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effectiveness of the maneuvers through the repeata- 
bility of the predicted parameters. Second, the com- 
parison of the average standard errors obtained from 
the computer programs with the standard error for 
the ensemble gives us a sense of how realistic the 
program values are. 

The 45 repeat maneuvers were analyzed using 
both the linear regression and the maximum like- 
lihood methods, and the resulting mean values of 
coefficients are compiled in table IV. Accompany- 
ing each parameter estimate is the average standard 
error associated with it as given by the particular 
estimation algorithm, the standard error for the en- 
semble of values, and the ratio of the two. The ratio 
given for each parameter is the factor by which the 
program-estimated standard error for that parame- 
ter should be multiplied to produce a realistic, ex- 
pected standard error. The ratios from the maximum 
likelihood results are seen to be greater than those 
from the linear regression. Two reasons are postu- 
lated for this difference: First, the standard error for 
each ML estimate is given by a numerical approxi- 
mation to the Cramer-Rao bound. Theoretically, the 
Cramer-Rao bound for an estimated parameter rep- 
resents the best achievable accuracy for that param- 
eter. The approximation to the Cramer-Rao bound 
used in the algorithm that produced the ML param- 
eter estimates presented in this paper is given by the 
diagonal elements of the inverse of an approximation 
to the Fisher information matrix. In general, for ac- 
tual flight data, these theoretical values are smaller 
than the estimated standard error from an SR al- 
gorithm. This gives, in general, smaller numbers in 
the denominators for the maximum likelihood ratios 
as compared with the linear regression ratios. Sec- 
ond, the dispersion of values for the maximum like- 
lihood estimates from these data was large for some 
parameters ( C Zq and C Zs , in particular). This type 
of dispersion is usually cfue to multiple correlations 
among variables in the data and unmodeled dynam- 
ics (such as those giving rise to the structure of the 
angle- of- attack residuals discussed in the section on 
data compatibility). 

All further analysis of flight data was principally 
carried out using the linear regression method be- 
cause of the consistency of results obtained through 
it, whereas some comparison estimates were made 
using the maximum likelihood method. 

Effects of Artificial Ice Shape 

Following the 45 repeat maneuvers completed 
with the aircraft in each of the 4 configurations, the 
pilot executed several additional longitudinal dou- 
blets at varying initial airspeeds (ranging from 62 to 


120 knots). These data were analyzed with the lin- 
ear regression and maximum likelihood algorithms to 
study the variations of parameters with initial flight 
condition and between different configurations. The 
results obtained from the stepwise regression routine 
are presented graphically in figures 5 through 8. In 
no case did the stepwise regression select a nonlinear 
term as a component of the aerodynamic model. 

Each parameter value is plotted with ±2a error 
bars, based on the standard error as given by the 
computer program. For the two configurations in 
which the flaps were held at 0°, there are statisti- 
cally significant differences between iced and uniced 
flights. These differences generally decrease with in- 
creasing velocity. When the flaps are extended and 
held at 10°, however, the distinction between iced 
and uniced becomes less dramatic. Disturbances in 
the flow coming off the wings caused by the extended 
flaps appear to negate a large percentage of the ef- 
fects of artificial ice attached to the leading edge of 
the tail section. The pitching-moment derivatives in- 
dicate that the artificial ice shape induces a loss of 
elevator effectiveness (C™ ), a decrease in static sta- 

bility {C ma ), and lower pitch damping (C m? ). 

Derivative values obtained from selected maneu- 
vers using the maximum likelihood routine are pre- 
sented in figure 9. Because of the convergence 
problems mentioned earlier, the maximum likelihood 
algorithm was run only on several selected maneu- 
vers. The resulting values were usually in agreement 
with the regression-determined parameters with re- 
spect to trends, but they differed significantly in ac- 
tual values estimated. No reason was found for those 
differences. 

Error Analysis 

A fundamental issue throughout this research is 
that of distinguishing between the effects of ice and 
the inherent variability of predicted parameters. In 
addition to the error analysis performed on the val- 
ues obtained from the 45 repeat maneuvers, the C m 

& e 

estimates for the 0° flap condition were considered 
(fig. 10) for an additional analysis. This set of points 
was chosen for this portion of the analysis because 
each of the two groups of values (iced and uniced) 
lend themselves to linear approximations. Also, de- 
termining a prediction interval (denoted by the terms 
in brackets) in this particular case serves to quantita- 
tively enforce what is qualitatively a very noticeable 
difference in parameters for the iced/uniced cases. 
After approximating the two groups of points by 
straight lines, prediction intervals about these lines 


were calculated as follows: 


[s/o(zo)]+ [< a ./2,»- 2 ] xs[l + n 1 +(*o -*) i / 5 'xx] 1/! 

( 3 ) 


where 


y o(*o) 

linear approximation of 
Cm 6 corresponding to some 
x$ value 

a * 

0.05 to achieve a 95-percent 
prediction interval 

n 

number of data points 

ta*/2,n-2 

student ^-distribution 

s 

standard error (can be one 
of three substitutions listed 
below): 

s pa (y) 

average of standard errors 
obtained from computer 


program 

*e(v) 

ensemble-average standard 
error for values of CV , 

® e 


[Uvi-y) 2 /(n-2)] l/2 

s r{v) 

s Pa(v) multiplied by ratio of 
ensemble to program average 
obtained from table IV 

Sxx 

sum of squares of residuals, 


N 

£(**-*) 2 

i=l 


The resulting interval is one within which we 
can expect 95 percent of all predicted elevator- 
effectiveness values to fall. Figure 10 presents these 


prediction intervals and gives an illustration of both 
the areas in which icing effects are discernible as well 
as the speed range in which we cannot attribute any 
change in Cm. to the artificial ice shape. A similar 
analysis can be done on other parameters by fitting 
them in a piecewise linear fashion. 

Conclusions 

Based on the analysis and results presented, the 
following four basic conclusions are presented: 

1. Any future research of this type will require 
determination of the cause of error found in angle-of- 
attack measurements. 

2. A multiplicative relationship between ensem- 
ble standard error and estimated standard error from 
the stepwise regression routine was identified with 
numerical values representing this ratio assigned to 
each stability and control derivative. 

3. A multiplicative relationship between ensem- 
ble standard errors and those obtained as an approx- 
imation to the Cramer-Rao bound as given by the 
maximum likelihood program was determined as a 
quantitative ratio for each parameter estimated. 

4. The effects of the artificial ice shape attached 
to the tail section have been shown to be measurable 
as changes in the longitudinal stability and control 
derivatives in all the flight test conditions addressed 
here with the largest effects being seen in those 
derivatives most directly determined by the flow on 
the tail, i.e., static stability ( C mQ ), pitch damping 
(Cm,), an d elevator effectiveness ( Cm 6e )• 

NASA Langley Research Center 
Hampton, VA 23665-5225 
January 25, 1989 
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Appendix 

Postulated Models 

The six-degrees-of-freedom equations of motion 
for the aircraft are 


Here, equations (4), (6), (8), and (10) describe the 
longitudinal motion and equations (5), (7), (9), and 
(11) describe the lateral motion of the aircraft. The 
angles of attack and sideslip are, respectively, given 
by 


u — —qw + rv — g sin 0 + ~~—C x 

2m 


(4) 


a = tan l (w/u) (3 = sin 1 (v/V) 


v — -ru + pw + g cos 0 sin </> + ^ — C Y 

2m 


w = — pv + qu + g cos 0 cos <f) + ~~^~C Z 

2m 

^•’ r ( k TT 1 -) + Tr {r '- r 




. pV 2 Sb„ 
V) + -7^—C n 


along with the kinematic relations 

0 = q cos <t> - r sin <f> 


(5) 

( 6 ) 

(7) 

( 8 ) 
(9) 


( 10 ) 


The aerodynamic forces and moments are repre- 
sented by the coefficients C a (where a = X, Y, Z, /, m, 
or n). It is postulated that these coefficients can 
be written as a Taylor-series polynomial expansion 
about an equilibrium trim condition as 


t 

C a = ©a,0 4" ^ ®a,j ^ x j (12) 

j=l 

where 0 a? o and 0 a j are the unknown parameters. 
For a = X, Z, or m, the independent variables Xj 
are a, qc/ 2V, fi e , and their combinations; and for 
a — y, /, or n, they are /3, p6/2V, r&/2V, <§ a , <5 r , and 
their combinations. Furthermore, 


and 


(j> — p 4* (g sin 0 + r cos 0) tan 6 (11) 


A xj = Xj(t) — Xj 


(* = 0 ) 
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Table I. Physical Characteristics of Research Aircraft 


Mass, kg 

Inertia: 

I x , kg-m 2 

Iy i kg-m 2 

I z , kg-m 2 

Ixz, kg-m 2 

Wing: 

Area, m 2 

Aspect ratio 

Span, m 

Mean geometric chord, m .... 
Airfoil section (17-percent thickness) 

Horizontal tail: 

Area, m 2 

Airfoil section (inverted) 


Low 

High 

4150 

4600 

21300 

21790 

30 000 

31030 

44 990 

48 640 

1430 

1500 


39.02 


10.06 


19.81 


1.98 


“deHavilland High Lift” 


9.1 

NACA 63A213 


Table II. Estimated Bias and Scale-Factor Errors Present in Longitudinal-Data 
Compatibility Maneuvers for Flight 88-5 


Estimated 

parameter 

Longitudinal compatibility maneuvers for — 

Run 33 

Run 34 

Run 35 

b a , rad 

0.0038 ± 0.0006 

0.0092 ± 0.0007 

-0.0047 ± 0.0007 

b q , rad/sec 

0.00057 ± 0.00003 

0.000996 dh 0.2E-09 

0.000645 ± 0.3E-07 

b$, rad 

-0.0175 ± 0.0005 

-0.0173 ± 0.0003 

-0.0207 ± 0.0004 

b ax , 9 units 

0.123 ± 0.003 

0.089 ± 0.030 

0.002 ± 0.002 

b az , 9 units 

0.014 ± 0.001 

0.004 ± 0.002 

-0.056 ± 0.001 

^ a 

0.022 ± 0.005 

-0.054 ± 0.008 

0.128 ± 0.010 


Table III. Estimated Bias and Scale-Factor Errors Present in Longitudinal Doublets 
and Large- Amplitude Maneuvers for Flight 88-5 


Estimated 

parameter 

Average of 
45 repeat 
maneuvers 

Large-amplitude maneuvers for — 

Run 36 

Run 37 

Run 38 

6 a , rad 

b q , rad/sec .... 

6$, rad 

b ax , g units .... 
b az , 9 units .... 
^ a 




-0.00381 ± 0.00232 
0.00129 ± 0.000151 
-0.0131 ± 0.00253 
-0.0144 ± 0.0163 
-0.0875 ± 0.00717 
0.1132 ± 0.0227 
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S(@) 

(c) 


0.0169 

.386 

.0122 

.0191 

.4314 

.0131 


s e (®)/s(Q) 


2.9 
2.7 
2.3 

.8 

1.5 

1.9 


(b) Maximum likelihood 


Standard errors 


Parameter 



“Mean value of ©. 
b SE(&)- ensemble standard error. 
c s(0): average program standard error. 



^(0)/s(0) 
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Figure 3. Control surface input and aircraft response for typical 2-1 elevator doublet. 


.2 

Angle of 
attack, .1 

rad 
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+ Measured value 
- Calculated response 

.01 

Angle-of- 
attack 

residuals, 0 

rad 

-.01 

0 5 10 15 20 25 30 

Time, sec 

Figure 4. Results from compatibility check on angle of attack. 
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Initial velocity, knots 

Figure 6. Pitching-moment derivatives and ±2cr error bounds from modified stepwise regression for flaps at 0 
















Initial velocity, knots 

Figure 9. Comparison of modified stepwise regression results with maximum likelihood results. Clean aircraft; 
flaps at 0°. 
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Figure 10. Confidence intervals of 95 percent for elevator-effectiveness derivatives for iced and uniced aircraft. 
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